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interpret the data. 
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1 Introduction 



Neutrino oscillation experiments have produced "smoking guns" for non-zero 
neutrino mass in the solar neutrino deficit , in the excess of p(T>e, e'^)n re- 
actions in the LSND experiment |^ , and from the zenith-angle dependence of 
the electron-to-muon event ratio in the Super-Kamiokande (SK) data[^ (see 
also [|[ |[ H). The results of reactor neutrino experiments @] constrain the dis- 
appearance of 77g well enough to imply that the SK data are dominated by 
Vfj^ —>■ Vr (z7^ — !■ Vr) oscillatioHS, with only a minimal oscillation to electron-type 
neutrinos, since they show Fg's do not oscillate as readily as required by the SK 
data for z^^ Vf^ or 77^ Vf, oscillations to have a significant role. 

While the interpretation of the SK data in terms of neutrino oscillations is 
widely accepted, there were questions concerning the interpretation of the LSND 
data as evidence of — > oscillations. Some doubt existed that the standard 
solar model was accurate enough to support the conclusion that there was really 
a deficit of solar neutrinos. When the results of all solar neutrino experiments 
are considered, there is no scenario in which these data are compatible with the 
standard solar model unless the flux of i/g from the sun oscillates partially into 
other i^-flavors to which the experiments arc not sensitive. 

On 17 July 2001, however, the Sudbury Neutrino Observatory (SNO) collab- 
oration settled this issue. They published results from the direct measurement 
of the rate of the reaction (i(i^e, s~)pp from solar neutrinos The neutrino flux 
implied by these data was compared with that from neutral-current neutrino- 
electron elastic scattering data from SK. It was concluded that there is an active, 
non-electron-flavor neutrino component in the solar neutrino flux, and that the 
total flux of active neutrinos from the *B branch is in close agreement with the 
standard solar model of Bahcall and his co-workers . The standard solar model 
is thereby conflrmed, and the case for neutrino oscillations is now compelling. 

A remaining question is that of the LSND indication of 17^ oscillations. 
All attempts to incorporate these results in the same analysis with the solar 
neutrino and atmospheric neutrino experiments fail in any scenario involving 
only three neutrino flavors. 

Accepting that now both solar and atmospheric neutrino experiments give 
clear evidence for neutrino oscillations, there are only two conclusions to be 
drawn from the LSND data. Either the excess events from the reaction p(l7e, e~^)n 
in the LSND is due to phenomena other than 77^; — ^ Ve oscillation, or that there 
must exist a fourth generation of neutrinos "sterile" with respect to "normal" 
weak interactions Q. To insist on accepting one or the other of these options 
at the present time is to accept an unsubstantiated theoretical prejudice. This 
issue is still very much an open one. The MiniBooNE experiment will settle this 
issue; however, it will be a search for v^. ,not — > as in LSND. 

While an unambiguous interpretation of all of the neutrino oscillation experi- 
ments is not yet possible, it is abundantly clear that neutrinos exhibit properties 
not included in the standard model, namely mass and flavor mixing. Accord- 
ingly, sensitive searches for neutrinoless double-beta decay {Qv /3/3-decay) are 
more important than ever. Experiments with kilogram quantities of germa- 
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nium, isotopically enriched in ^^Ge, have thus far proven to be the most sen- 
sitive, specifically the Heidelberg-MoscowJlO| and IGEX||l^ experiments. The 
lower limits in half-life sensitivities are: 1.9 x 10^'°^ y[0 and 1.6 x 10^^ y||ll|- A 
new generation of experiments will be required to make significant improvements 
in sensitivity, as discussed later. 

Petcov and Smirnov|T^ showed that both MSW and vacuum oscillation so- 
lutions of the solar neutrino problem can be compatible with Ov /3/3-decay driven 
by an effective Majorana electron-neutrino mass in the range 0.1 to 1.0 eV. The 
interpretation of all the neutrino oscillation data together, as discussed later, 
implies a range that could be between 5 and 10 times lower. The exploration 
of such a range will require next-generation experiments. Some that are being 
proposed are: CAMEO(l|, CU0RE114|, EX0||, GENIUS|6|, Majorana0, 
and MOOnIisII. 



2 Neutrinoless Double-Beta Decay 

Neutrinoless double-beta decay is a process by which two neutrons in a nu- 
cleus beta decay by exchanging a virtual Majorana neutrino, each emitting an 
electron. This violates lepton number conservation (AZ = 2). There are many 
reviews on the subject [|l^, |2l|, |2^ . 

The decay rate for the process involving the exchange of a Majorana neutrino 
can be written as follows: 

X% = G°^(£;o, Z){m^)'\M'}'^ - {9a/ 9v)Ml^r,\\ (1) 

In equation |l|, G^'^ is the two-body phase-space factor including coupling 
constants, M^^ and Mg,^ are the Fermi and Gamow- Teller nuclear matrix ele- 
ments respectively, and qa and gv are the axial-vector and vector relative weak 
coupling constants, respectively. The quantity (mi^) is the effective Majorana 
neutrino mass given by: 

2n 

{m,)^\Y^X^^{Utu?mul (2) 

k=l 

where A^^ is the CP eigenvalue associated with the k*^ neutrino mass eigenstate 
(±1 for CP conservation), J7;^ is the (Z,fc) matrix element of the transforma- 
tion between flavor eigenstates and mass eigenstates \vk) for left handed 
neutrinos, 

\vi)^Y.^tkW). (3) 

and m/c is the mass of the k*^ neutrino mass eigenstate. 

The effective Majorana neutrino mass, (mi,), is directly derivable from the 
measured half-life of the decay as follows: 

(m,) = vi,{FnT^,J^)-^/\V, (4) 
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where Fn = G^^IM^" - [qa/ gv)M^j,\^ , and me is the electron mass. This 
quantity derives from nuclear structure calculations and is model dependent. 

The most sensitive experiments thus far utilize germanium detectors isotopi- 
cally enriched in ''^Ge from 7.78% abundance to ~86%. Germanium detector 
experiments were started with natural abundance detectors by Fiorini, et al., in 
Milan|Q , evolving over the years to the first experiments with small isotopically 
enriched detectors, and finally to the two present multi-kilogram isotopically- 
enriched ^^Ge experiments: Heidelberg-Moscow [p^ and IGEX|pl|. Reference 
[ p!o| has about four times the exposure as reference with limits of similar 
magnitude. This strongly implies that experiments with on the order of 100 
moles of ^^Ge have reached the point of diminishing returns. 

Where should the field of double-beta decay go from here? Suppose the 
observed neutrino oscillations of atmospheric neutrinos and solar neutrinos are 
considered. What probable range of {va.^) is implied? Would it be large enough 
to allow the direct observation of Qv /3/3-decay? If so, what technique would be 
the best for a possible discovery experiment? 

3 Theoretical Motivation and Probable Neutrino 
Scenarios 

The SK data imply maximal mixing of with v-r with Siri^r^ ~ 3 x 10~^ eV^. 
The solar neutrino data from SK and from SNO also imply that the small mixing 
angle solution to the solar neutrino problem is disfavored, so that 5m^ (solar) 
~ (10~^ — 10"^) eV^. Based on these interpretations, one probable scenario for 
the neutrino mixing matrix has the following approximate form: 

Ve \ ( 1/V2 l/\/2 \ / 1^1 \ 

= -l/\/2 1/V2 1/V2 ^2 (5) 

Vr ) \ 1/V2 -1/V2 1/V2 / \ / 

The neutrino masses can be arranged in two hierarchical patterns in which 
(5m|i ~ (5m§2 3 x 10"^ eV^, and Sml-^ - (IQ-^ - IQ-*) eV^. With avail- 
able data, it is not possible to determine which hierarchy, 7713 > 7711(7712), or 
'771(7772) > 7713, is the correct one, nor is the absolute value of any of the mass 
eigenstates known. 

The consideration of reactor neutrino and atmospheric neutrino data to- 
gether strongly implies that the atmospheric neutrino oscillations are very dom- 
inantly t'^ Vt(v^ — > v^), which implies, as seen from equation 5, that is 
very dominantly a mixture of vi and V2- In this case there will be one relative 
CP phase, e, and equation |^ reduces to the approximate form: 

(m^> = i(mi +em2), (6) 

where the large mixing angle solution of the solar neutrino problem implies 

(7772 - 777,2) ^ (^0-5 _ ^0-4) eV^. (7) 
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Assuming that e ~ +1, and that neutrinos are Majorana particles, then it is 
very probable that Ov /3/3-decay is driven by an effective electron neutrino mass 
between 0.01 eV and the present bound from ^^Ge experiments. Consider the 
relation 

0, {ln2)Nt 

-'1/2- > (P) 

where A'' is the number of parent nuclei, t is the counting time, and c is the 
maximum number of counts that can be attributed to Ov /3/3-decay. To improve 
the sensitivity to (m^) by a factor of 10^ from the present ~1 eV to 0.01 eV, 
one must increase the quantity Nt/c by a factor of 10^. The quantity N can 
feasibly be increased by a factor of ~ 10^ over present experiments, implying 
t/c must also be improved by a similar amount. Since the present counting 
times are probably about a factor of 5 less than a practical counting time, the 
background should be reduced by at least a factor of 20 below present levels. A 
further reduction is probably feasible and should be pursued. 



4 The Majorana '°Ge Oiy /3/3-Decay Experiment 

The Majorana experiment is proposed for a US deep underground laboratory, 
and requires very little R&D. It stands on the technical shoulders of the IGEX 
experiment and other previous successful double-beta decay and low-background 
experiments. Furthermore, new segmented Ge detector technology has recently 
become commercially available, while Pacific Northwest National Laboratory 
(PNNL)/University of South Carolina (USC) researchers have developed new 
pulse-shape discrimination techniques. 

The IGEX experiment terminated with 117 mole-years of data with an av- 
erage background of 0.06 c/keV/kg/y in the later data sets. The resulting half 
life is 1.6 X 10^^ y and its implied bounds on (mi,) are shown in Table 1. 

The Majorana experiment represents a great increase in Ge mass over IGEX 
with new segmented Ge detectors and the newest electronic systems for pulse- 
shape discrimination. It is conceived to be 500 kg of Ge detectors, fabricated 
from Ge isotopically enriched to 86% in ''^Ge. The typical detector size will 
be approximately 2.4 kg, requiring about 200 detectors. Below is presented 
the projected sensitivity for such an array, shown in an artist's conception in 
Figure |^. 

The array will have a fiducial mass of 500 kg, containing A^ = 3.43 x 10^^ 
atoms of ''^Ge; a counting time of 10 years, and an energy resolution of 3.0 keV 
FWHM. The starting background assumed is that achieved by IGEX, prior 
to pulse-shape discrimination, which was 6o = 0.2 counts/keV/kg/y at the 
Ov /3/3-decay energy of 2039 keV. This background was achieved at 4000 mwe 
and was completely accounted for as the decay of cosmogenic isotopes in the 
germanium||. The decay of ^^Ge (270.8-day half-hfe) and ™Co (5.7-year 
half-life) over the 10-year counting period, combined with estimates of reduced 
above-ground cosmogenic activation and additional decay during detector con- 
struction, reduces this background by an overall factor of about twenty, com- 
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pared to the IGEX background mentioned. This reduces the expected back- 
ground rate to 0.01 c/keV/kg/y. 

The optimum energy window for a simple Poisson analysis is 2.8(7 of the 
Gaussian peak (83.8%), or 3.57 keV. The total number of background counts in 
this window around 2039 keV over a 10-year period would then be 

b-AE-M-t= (0.01 c/keV/kg/y)(3.57 keV)(5000 kg x y) = 178 (9) 

prior to data cuts. 

Above-ground experiments by the PNNL/USC collaborators demonstrate 
that new pulse-shape discrimination methods have an acceptance fraction for 
internal cosmogenic backgrounds of no more than 0.265. Monte Carlo simula- 
tions demonstrate detector segmentation will allow an additional cut with an 
acceptance fraction for internal cosmogenic backgrounds of no more than 0.138. 
An example of the effect of segmentation on internal cosmogenic background 
in the 2 MeV region-of-interest is shown in Figure Applying both of these 
orthogonal cuts is expected to reduce the background to 6.5 remaining events. 
If the background in the region is featureless in energy, a simple Poisson analysis 
yields a limit of 6 events as the number of /9/3-decay events consistent with 
the background (90% CL). 

The two cuts discussed do not have 100% efficiency for accepting 13(3- 
decay events. The pulse-shape discrimination has a measured single-site-event 
acceptance fraction of 0.802. The detector-segmentation cut has a calculated 
single-site-event acceptance fraction of 0.907. The total counting efficiency is 
then 0.727. Accordingly the sensitivity is projected as: 

T^^sens) . (0-727)(.693)(3.43 x 10^^) 10 ^ ^ ^ ^^^^ 

This analysis does not account for the fact that the close-packed granularity 
of ~200 crystals will allow further cuts against multiply-scattered background 
events. The background computed above is therefore conservative. The result- 
ing upper bounds on (m,y) are shown in Table |l| and, discarding an unfavored 
QRPA form factor, range from 0.02 eV to 0.07 eV. 

Table 1: Nuclear structure factor F/v and Majorana neutrino mass parameters 
(m^) for T°J^ = 1.6 x lO^^ y and for 4.2 x 10^^ y. 



Fj^,yT 1 Model [reference 


(m,) (eV) (m,) (eV) 


1.56 X 10-1^ SheU Model 
9.67 X 10-15 QRPA 
1.21 X 10"" QRPA 
1.12 X 10-" QRPA 
1.41 X 10-" SheU Model 


25 
26 
27 
28 
29 


0.32 0.020 
1.3 0.080 
0.37 0.023 
0.38 0.024 
1.08 0.067 
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Figure 1: An artist's conception of the Majorana array of Ge detectors using 
present cryostat technology. Other, more granular arrays will also be investi- 
gated. 
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Figure 2: Monte Carlo simulation of internal ^"Co background. The plot on the 
left shows the spectrum before and after a one-segment-only cut is applied. The 
plot on the right displays a histogram of the number of segments in which energy 
was deposited for internal ^°Co background events falling in a 2.0-2.1 MeV 
region-of-interest . 
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5 Recent Progress in Ge Detector Technology 



Majorana is not simply a volume expansion of IGEX. It must have superior 
background rejection. Because it has been conclusively shown that the lim- 
iting background in at least some previous experiments has been cosmogenic 
activation of the germanium itself, it is necessary to mitigate those background 
sources. Cosmogenic activity fortunately has certain factors which discrimi- 
nate it from the signal of interest. For example, while Qi' /3/?-decay would 
deposit 2 MeV between two electrons in a small, perhaps 1 mm'^ volume, in- 
ternal ^°Co decay deposits about 318 keV (endpoint) in beta energy near the 
decaying atom, while simultaneous 1173 keV and 1332 keV gammas can deposit 
energy elsewhere in the crystal, most probably both in more than one location, 
for a total energy capable of reaching the 2039 keV region-of-interest. A similar 
situation exists for internal ^^Ge decay. Thus deposition- location multiplicity 
distinguishes double-beta decay from the important long lived cosmogenics in 
germanium. Isotopes such as ^^Co, ^"^Co, ^^Co, and ^^Ge are produced at a rate 
of roughly 1 atom per day per kilogram on the earths surfacej2^. Only ^°Co 
and ^*Ge have both the energy and half- life to be of concern. 

To pursue the multiplicity parameter, two approaches are possible. First, the 
detector current pulse shape carries with it the record of energy deposition along 
the electric field lines in the crystal; crudely speaking, the radial dimension of 
cylindrical detectors. This information may be exploited through pulse-shape 
discrimination. Second, the electrical contacts of the detector may be divided 
to produce independent regions of charge collection. 

By segmenting the inner contact into two (axial) parts and the outer contact 
into 6 (azimuthal) parts, as shown in Figure ^ multiplicity data can be obtained. 
The Monte-Carlo simulation data set shown in Figure ^ is based on this con- 
figuration and shows that internal highly- multiple backgrounds like ^"Co can 
be strongly suppressed at 2038 keV. The internal ^"Co modeled in the figure is 
produced by cosmic-ray neutrons during the preparation of the detector, accu- 
mulating after the last crystal-growth step. Its elimination by segmentation and 
pulse shape discrimination is crucial. Beyond this simple segmentation cut, it 
may be possible to use the signals derived from segments seeing no net charge, 
adjacent to a segment seeing net charge, to locate a single-site deposition in 
the axial and azimuthal coordinates of the crystal or to distinguish a single-site 
deposition from a multiple one. This possible technique is not currently in- 
cluded in sensitivity calculations. Intrinsic n-type Ge detectors with segmented 
electrical contacts are now available commercially. After Monte Carlo studies 
and discussions with detector manufacturers, the general configuration shown 
in Figure |^ is believed to be optimal for balancing background reduction, cost 
and production efficiency. 

The development of the pulse-shape discrimination technique used in the sen- 
sitivity calculation arose from earlier work comparing experimental pulses with 
completely simulated ones. After repeated trials with fully simulated pulses, 
this early method was set aside due to the unpredictable effect of preamplifier 
noise levels and the unacceptable sensitivity of the resulting cuts to the detector 
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Twelve Section Segmented Ge Detector 



8 cm 




Top View 

0.9 cm 



8 cm 



Side View 

Figure 3: The geometry of the segmented detectors planned for the Majorana 
experiment. 



and charge-collection model. For example, the parameters needed to correctly 
simulate pulses include the space-charge density remaining after crystal deple- 
tion. This is most often poorly known and represents essentially a set of free 
parameters for the simulation. The operating voltage of the crystal, which may 
require variation during the crystal lifetime, is also critically important to the 
fidelity of simulated pulse formation. Thus, any system built upon the accuracy 
of predicted or modeled pulse-shapes was found less than optimum for cutting 
efficacy and fragile with respect to long-term maintenance. 

A much-superior pulse-shape discrimination method was constructed using 
a set of parameters calculated from each pulse. The parameters for a single 
pulse are compared to a distribution collected from a calibration source of in- 
terest. Those which conform to a large degree are kept; those which do not 
are discarded. The acceptance ratio is adjusted during pulse post-processing 
to optimize cut efficacy. The strength of this approach is that no set of simu- 
lated pulses need be maintained; no space charges need be known. It is only 
required that calibration data be obtained before and after any voltage change, 
a procedure which is required in any case to periodically check detector gain. 
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To summarize, new techniques depend entirely on experimental calibration and 
do not utilize pulse-shape libraries. The ability of these techniques to be eas- 
ily calibrated for individual detectors makes them practical for large detector 
arrays. 

This development was in large part due to a shift in thinking about digi- 
tization, earlier digitization had been done with an 8-bit digitizer sampling as 
frequently as 500 MHz. Theory held that by oversampling, the energy resolu- 
tion afforded by analog acquisition could be achieved, which is well-matched to 
14-bits. However, measurement of actual preamplifier bandwidth showed that 
there is no more than 25 MHz of information in the signal. Thus, emerging 
hardware for 40-MHz 12-bit acquisition, such as the Digital Gamma Finder 
(DGF) produced by X-Ray Instrumentation Associates (XIA)|Q, allows both 
a high-resolution energy value and a high-fidelity waveform to be acquired in 
a single, lower cost, and highly stable electronics chain. After acquiring a de- 
tector preamplifier pulse, all subsequent operations on the signal are performed 
digitally. The particular unit used in pulse-shape discrimination development is 
the 4-channel DFG-4C, developed and manufactured by XIA. 

The DGF-4C has 4 independent, 40 MHz, 12-bit analog-to-digital converters 
(ADCs). The ADCs are followed by First-in, First-out (FIFO) buffers capable 
of storing 1024 ADC values for a single event. In parallel with each FIFO is 
a programmable digital filter and trigger logic. The digital filter and trigger 
logic for each channel is combined into a single Field Programmable Gate Array 
(FPGA). Analog input data are continuously digitized and processed at 40 MHz. 

Experimental example pulses are shown in Figure ^. An example single-site 
event from the 1592 keV double-escape peak of the ^osxi 2615 keV line is shown 
as the top signal. The bottom signal is an example multi-site pulse from the 
full-energy peak of the ^^^Bi line at 1620 keV. 

The DGF-4C is a smart filter of incoming pulses. If, for example, a signal 
has a pulse-width incompatible with the usual collection time of 200-1000 ns, 
or is oscillatory (like microphonic noise), the trigger filter can be programmed 
to reject it. This feature can be used to allow the very low energy thresholds 
required in Dark Matter searches as well as eliminating a broad spectrum of very- 
low-rate artificial pulses from high voltage leaks and electromagnetic interference 
that can appear in the energy region of Ov /3/?-decay. 

Calibration for single-site event pulses is trivially accomplished by collecting 
pulses from thorium ore; the 2614.47-keV gamma ray from ^os^pj produces a 
largely single-site double-escape peak at 1592.47 keV. The PSD discriminator is 
then calibrated to the properties of the double-escape peak A slightly improved 
double-escape peak can be made from the ^^Al gamma ray of 2938.22-keV. 
The double-escape appears at 1916.22 keV, only about 120 keV away from the 
expected region of interest for Oi/ /3/3-decay. The obvious and direct use of pulse- 
shape discrimination and segmentation is the rejection of cosmogenic pulses in 
the germanium itself. However, the approach should be also effective on gamma 
rays from the shielding and structural materials. 

The effects as background of neutrons of both high energy (cosmic muon 
generated) and low energy (fission and {a,ii) from rock) are under consideration 
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Smoothed Current Pulse (SG 9 points, order 4) 
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Figure 4: Reconstructed current signal from a 1592-keV double-escape-peak 
event (top) and a 1620.6-keV full-energy gamma peak event (bottom). Hori- 
zontal scale spans 2500 ns. 



as well. The segmentation and granularity of the detectors will provide some 
protection from this lower-order background. These neutrons could also produce 
other unwanted activities. For instance, and "'^''C can be produced in nitro- 
gen from high and low energy neutrons, respectively. Fortunately, Majorana 
detectors will not be surrounded by nitrogen at high density. 

In conclusion, the Majorana project has been designed in a compact, mod- 
ular way such that it can be built and operated with high confidence in the 
approach and the technology. The initial years of construction will allow al- 
ternate cooling methods to be employed if they have an advantage and should 
they be shown to overcome long-term concerns due to surface contamination, 
muon-induced ions, and diffusion. The technology supporting Majorana signal- 
processing-based background rejection is ready and shows promise of future 
improvements. The technology of radiopure copper production has improved 
steadily since the last IGEX system, in which copper support materials played 
no measurable radiological role. Finally, those who enrich germanium have also 
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expressed their enthusiasm and support for the project. Thus, the Majorana 
project is poised ready to begin operations toward determining the effective 
Majorana mass of the electron neutrino. 
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